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Abstract 

We studied a constraint on the primordial helium abundance Y p from current 
and future observations of CMB. Using the currently available data from WMAP, 
ACBAR, CBI and BOOMERANG, we obtained the constraint as Y p = 0.25±g;J? at 



C.L. We also provide a forecast for the Planck experiment using the Markov 
chain Monte Carlo approach. In addition to forecasting the constraint on Y p , we 
investigate how assumptions for Y p affect constraints on the other cosmological pa- 
rameters. 



1 Introduction 



Current cosmological observations are very precise to give us a lot of information on the 
evolution and present state of the universe. We usually extract the information by con- 
straining cosmological parameters such as energy densities of baryon, dark matter and 
dark energy, the Hubble constant, reionization optical depth, spectral index of primordial 
fluctuation and so on. Among them, in this paper, we focus on the primordial helium 
abundance Y p , which has been of great interest in cosmology. One of the reasons why the 
primordial helium abundance has been considered to be interesting and important is that, 
in the context of the standard big bang nucleosynthesis (SBBN), we can know the baryon 
density once Y p is determined from observations. However, it has been discussed that a 
significant systematic error dominates when one infers the value of Y p from measurements 
in low-metallicity extragalactic HII region [1-5]. Furthermore, there have been some dis- 
cussions that there may be a large uncertainty in the neutron lifetime [6-8] , which results 
in uncertainties in the predictions for the abundances of light elements. In this respect, the 
study of other independent measurements of the helium abundance would be interesting. 

Recent precise measurements of cosmic microwave background (CMB) such as WMAP 
can now enable us to constrain cosmological parameters with great accuracies. However, 
the helium abundance has not been discussed much in the literature when one study of 
cosmological constraints from CMB since Y p has been considered to have little effect on 
CMB power spectrum. In most of analyses, Y p is fixed to be 0.24 which is probably 
motivated from a somewhat old value of the observed primordial helium abundance of 



Y p = 0.238 [9JO But, in fact, there have been some works which discuss the effects of 
Y p on CMB and some constraints were given [10-12]. Since the helium abundance affects 
the recombination history, the CMB power spectrum can be affected mainly through the 
diffusion damping. Although the constraint is not so severe, it is important to notice 
that they are obtained independently from BBN, which can be used to cross-check our 
understanding of the helium abundance. Furthermore, the value of Y p at the time of BBN 



may be different from that at late time when CMB observations can probd^ 2 !. After we 
studied the constraint on Y p in [12], the data from WMAP has been updated [14-18]. In 
addition, the data at higher multipoles where the effects of Y p become significant have 
been updated by ACBAR [19] and CBI [20]. Thus it is a good time to investigate the 
constraint on the helium abundance using these CMB data, which is one of the aims of 
this paper. 

Furthermore, we expect a more precise measurement of CMB from the future Planck 
satellite [21], which can give us a much better constraint on Y p . In fact, future constraint 
on Y p has already been studied using the Fisher matrix formalism [10,12]. Although this 
method is fast and usually adopted to predict the precision of the future measurements of 



#1 Recent observations give e.g. Y p = 0.249 ± 0.009 [1,2]. 

# 2 For example, in a scenario to solve so-called "Lithium problem" with Q-balls, the decay of Q-balls 
produce extra baryon after BBN has completed [13]. In a model of this kind, the value of Y p can vary at 
different epochs. 
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cosmological parameters, it can give some inaccurate predictions in some cases, for exam- 
ple, when the likelihoods does not respect a Gaussian form. In addition, the Fisher matrix 
formalism predicts only the uncertainty for the parameter estimation since it just con- 
cerns with the derivatives with respect to parameters around the fiducial values. However, 
since some parameters are correlated in general, fixing the values of some parameters can 
bias the estimation of other parameters. Namely, priors we assume on some parameters 
can cause the estimated central values to deviate from the input fiducial values, but such 
effects cannot be quantified by the Fisher matrix approach. Thus, in this paper, we use 
the Markov chain Monte Carlo (MCMC) approach to extract reliable future constraints 
on Y p and other cosmological parameters. In particular, when we forecast the sensitivity 
for other cosmological parameters, we assume some different priors on Y p and investi- 
gate to what extent the information of Y p is important to determine other cosmological 
parameters. 

The organization of this paper is as follows. In the next section, the effects of Y p on 
CMB power spectrum are briefly discussed. Then in section [31 we study the constraint 
on Y p using CMB data currently available including ACBAR, BOOMERANG and CBI as 
well as WMAP5. In section HI we investigate a constraint on Y p and other cosmological 
parameters from the future Planck experiment. In addition, we also study how the prior 
on Y p affects the determination of other cosmological parameters. A brief discussion on 
the significance of the uncertainties of the recombination theory in deriving cosmological 
constraint is given too. The final section is devoted to the summary of this paper. 

2 Effects of Y p on CMB 

Here we briefly discuss the effects of the primordial helium abundance Y p on CMB power 
spectrum where Y p = An^e/ (^h + 4riH C ) with tih and nn c being the number density of 
hydrogen and helium-4 respectively. As has been discussed in [10,12], the value of Y p 
can affect the recombination history, which changes the structure of acoustic peaks. The 
effects of Y p on acoustic peaks mainly come from the diffusion damping which suppresses 
the power on small scales and the shift of the position of acoustic peaks due to the change 
of the recombination epoch. Before recombination, the number density of electron n e 
can be given by n e = n&(l — Y p ), where ni, is the baryon number density. Thus the 
increase of Y p indicates the decrease of the number of electrons. When the number of 
electrons is reduced, the mean free path of the Compton scattering becomes larger, which 
means that fluctuations on larger scales can be more affected by the diffusive mixing and 
rescattering. Thus the damping scale below which fluctuation of photon is exponentially 
suppressed becomes larger. Furthermore, due to the change of the number density of 
electrons, the epoch of recombination is also affected even though its effect is not so 
significant. This effect shifts the position of acoustic peaks slightly. In Fig. [TJ we show 
the CMB TT spectrum with several values of Y p . For reference, we also plot the current 
data (left panel) and the expected data from the future Planck experiment (right panel). 
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As mentioned above, by increasing the value of Y p , the power on small scales is damped 
more. In addition, it is noticeable that the position of acoustic peaks also shifts slightly. 




Figure 1: CMB TT power spectra for several values of Y p . In this figure, we take Y p = 0.1 
(blue dotted line), 0.24 (red solid line) and 0.4 (green dashed line). Other cosmological 
parameters are assumed to be the mean value of WMAP5 for a power-law ACDM model. 
For reference, in the left panel, data from WMAP5, ACBAR, BOOMERANG and CBI 
are also depicted. In the right panel, the expected data from the Planck experiment are 
also shown. 



To characterize the effects of the change in Y p and other cosmological parameters on 
CMB TT power spectrum Cf T , we consider some useful quantities [22]. First of all, to 
see how cosmological parameters affect the position of acoustic peaks, we investigated 
the response of the position of the first peak by the change of the parameters, which we 
denote All. In addition, to see the effects of the diffusion damping and some other effects 
by cosmological parameters, we study the height of the first peak relative to that at / = 10 
and the height of the second peak (and higher peaks up to 5th peak) relative to the first 
peak, which are denoted as Hi, H 2 , H 3 , if 4 and if 5, respectively. For clarity, we give the 
definitions of these quantities. The definition of Hi is 

_ / Ar(f = fo y 
Hi = {at(i = io)J ' (1) 

and the height of the i-th peak relative to the first peak is defined as 

(AT(l = U) 



where (AT(Z)) 2 = l(l + l)Cf T /2ir. We varied cosmological parameters including Y p around 
a fiducial model and obtained partial derivatives by fitting linearly around the fiducial 
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value. For the fiducial cosmological values, we assumed the mean values of the WMAP 
5-yr result (WMAP5) for a power-law ACDM model. Regarding Y p , we take Y p = 0.248 
which corresponds to the value obtained in the SBBN for the WMAP5 baryon density. 
The resulting derivatives are: 

Aujh „ Au m An s AY,, Ah , , 

A/i = 15.6 - 27.0 + 36.0 + 0.94-/ - 44.5—, (3) 

Uf, u m n s Y p h 

Aujh Au m An s Ah , , 

A#i = 2.87 - 3.13 - + 16.7 - - 2.30— , 4 

u b u m n s h 

Auji, Auj m An* AY n 
AH 2 = -0.290 + 0.023 + 0.396 - 0.013 (5) 

u b u m n s Y p 

Aujh Auj m An s AY V 
AH 3 = -0.177 + 0.206 + 0.514 - 0.028—/, (6) 

Uf, u m n s Yp 

Aujh Au m An* AY n , , 

AH 4 = -0.102 + 0.082 + 0.317 - 0.025—/, (7) 

Uf, u m n s Yp 

Aujh Auj m An, AY„ , x 

AH, = -0.040 + 0.084 + 0.236 - 0.023- p 



LUh tu m n s Y % 



p 



where cuj, and u m are energy densities of baryon and matter, n s is the scalar spectral 
index of primordial fluctuation, h is the Hubble constant in units of 100 km s _1 Mpc~ 1 . In 
the formula for Hi, we do not show the dependence on AY P /Y P since its effect on Hi is 
very small compared to that of the other parameters. As seen from the negative signs of 
AHi/AYp for % = 2-5, the diffusion damping becomes more efficient as Y p increases. We 
can also see the correlation of Y p with other cosmological parameters, which can be useful 
when we interpret the results, in particular, for a Planck forecast. 



3 Current constraint on Y p 

Now we discuss the constraint on Y p from current cosmological observations. For this pur- 
pose, we make use of the CMB data from WMAP5 [14-18], ACBAR [19], BOOMERANG 
[23-25] and CBI [20]. To investigate the constraint, we performed a MCMC analysis by 
using a modified version of cosmomc code [26]. We sampled in an 8 dimensional parameter 
space with (a>&, uj c , t, 6 s , n s , A s , Y p , Asz) where uj c is the energy density of dark matter, r 
is the optical depth of reionization, 9 S is acoustic peak scale [27], A s is the amplitude of 
primordial curvature fluctuation at the pivot scale k Q = 0.05 Mpc -1 and Asz is the ampli- 
tude of thermal Sunyaev-Zel'dovich (SZ) effect which is normalized to the Cf z template 
from [28@. In this paper, we consider a flat universe and assume a cosmological constant 
as dark energy. We also assume no running for primordial scalar fluctuation and no tensor 
mode. When we report our results in the following, we also use other customarily used 



# 3 However, the SZ effect may be so farge at very high multipoles that this template may not be 
appropriate to adopt. Hence, we conservatively do not use the ACBAR and CBI data with I > 2100. 
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prior ranges 


parameters 


current data 


Planck 


UJb 


0.005 


— > 0.1 


0.005 — > 0.1 


u c 


0.01 - 


* 0.99 


0.01 — > 0.99 




0.5 - 


-> 10 


0.5 -> 10 


T 


0.01 - 


-> 0.8 


0.01 -> 0.8 


n s 


0.5 - 


* 1.5 


0.5 -> 1.5 


lntlO 10 ^) 


2.7- 


-> 4 


2.7^4 


Y P 


(o- 


-1) 


(o-i) 


A sz 


- 


+ 2 




F H 






(0-2) 


bHe 






(0-1.5) 



Table 1: The prior ranges for the parameters used in the analysis. Priors shown in the 
1st and 2nd columns are adopted for current and expected Planck data, respectively. Note 
that we adopt the prior range for Y p shown above only in the cases with Y p being treated 
as a free parameter whereas Y p is a derived parameter in the case where we assume the 
SBBN relation. For the analysis with the current data, we also vary the amplitude of 
the Sunyaev-Zel'dovich effect Asz, which is omitted in the Planck data analysis. We also 
include two additional parameters Fh and bue which represent uncertainties in the theory 
of recombination (see section H] for more details). 



cosmological parameters such as the Hubble constant Hq = 100k km s Mpc~ and energy 
density of matter Q m = (ub + uj c )/h 2 . In performing a MCMC analysis, we impose top-hat 
priors on the primary parameters given above, which are summarized in Table [TJ 

Now we discuss the constraints on Y p when WMAP5 alone is used and when the 
data of ACBAR, BOOMERANG and CBI are used in addition. In Fig. ^ we show 
one-dimensional marginalized distributions of Y p for these two cases and, in the 1st and 
2nd columns of Table [2], we list the parameter estimations for Y p and other cosmological 
parameters. When we use WMAP5 data alone, the constraint is given as Y p < 0.44 at 
95% CL. where we only give the upper bound since the likelihood has a sizable value 
at Y p = 0. On the other hand, when all data are combined, the constraint is Y p — 
0.25^;JJ^J;J^ at 68% (95%) CL. For the analysis using WMAP5 alone, the limit we 
obtained is consistent with the result given in [15]. We see that by including the data 
from ACBAR, BOOMERANG and CBI the likelihood distribution has a well-defined 
peak which is close to Gaussian. It may also be interesting to notice that when the data 
from ACBAR and BOOMERANG and CBI are combined, the mean value becomes as 
Y p = 0.25, which is very close to the value obtained from HII region observations although 
the uncertainty is still large. With the help of this data set, we may begin to see the 
concordance with regards the Y p measurement from CMB and that from the HII regions. 
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Figure 2: One-dimensional marginalized distribution of Y p for the cases with WMAP5 
alone (red solid line) and all CMB data combined (green dashed line). 



Next we discuss the effects of the prior of Y p on the determination of other cosmological 
parameters. For this purpose, we repeated a MCMC analysis fixing the value of the helium 
abundance to Y p = 0.24 as in usual analyses. We use all the CMB data (i.e., WMAP5, 
ACBAR, BOOMERANG and CBI) here. In Table [2], in the last column, the constraints on 
cosmological parameters for the case with fixing Y p = 0.24 are shown. When we compare 
the constraints for the cases with and without fixing Y p , the central values as well as the 
errors at 68% C.L. are almost unchanged. Thus we can conclude that the usual practice 
of fixing of Y p = 0.24 scarcely affects the constraints on other cosmological parameters 
with the current precision of CMB data. However, since we can expect more precise 
measurements of CMB in the near future, the prior on Y p may become important and can 
affect the constraints on other cosmological parameters. We study this issue in the next 
section. 

4 Future constraint from Planck 

In this section, we forecast a constraint for the Planck experiment [21] focusing on the 
constraint on Y p itself and how the prior on Y p affects the constraints on other cosmological 
parameters. In fact, constraints from the future Planck experiment from this viewpoint 
have already been discussed in Refs. [10, 12] by using the Fisher matrix analysis. As 
mentioned in the introduction, when the likelihood of cosmological parameters can be 
approximated by a multivariate Gaussian function, the Fisher matrix analysis can give a 
reliable prediction. However, in practice, the likelihood function deviates from the Gaus- 
sian form. Furthermore, since the Fisher matrix analysis can only predict the uncertainty 
for a fiducial value, it cannot extract a bias effect (i.e. the estimated central value deviates 
from the fiducial value) which is caused by assuming priors on parameters and possible 
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parameters 


WMAP5 alone 
(Y p free) 


CMB all 

free) 


CMB all 
(Y p = 0.24) 


UJb 

u c 

T 

n s 

^(lO 10 ^) 

Y P 
A S z 


— U.UUb 

0.0881*™ 

3.06 ±0.06 
< 0.44(95%) 

i 1+0.9 
i - i -0.3 


n n99Q 4- n nnn^ 

n 1 1 0+0.006 
u.i±o_ 007 

1.043 ±0.004 

0-0871°;™ 

0.9671°;™ 

3 08+ - 04 
o.uc>_ 005 

25+ - 10 

u - zd -0.07 
i i +0.9 
i - i -0.3 


n229 +u - 0uUti 

U.U^Z3_ Q005 
n i 1 9+0.005 
u - llz -0 .006 

1.043 ±0.003 
0.0871°;™ 
0.9661°;™ 
3.07 ±0.04 

1 Q + 1.0 


Ho 


0.25 ±0.03 
72.3±1:1 


0.27 ±0.03 

>j-i y+2.2 
' ± -' -2.6 


u - zo -0.03 
71 y+2.3 
' L - ' -2.5 



Table 2: Mean values and 68% errors from current observations of CMB for the cases with 
WMAP5 alone and all data combined. (Regarding Y p , an upper bound at 95% C.L. is 
given for the case with WMAP5 alone.) In the last column, the value of Y p is fixed as 
Y = 0.24. 



correlations among parameters. Thus it may be better to make a more reliable prediction 
by using a MCMC method. For this purpose, we follow the approach of Ref. [29]. 

Here we briefly explain the method of Ref. [29] . Observed anisotropics can be expanded 
in spherical harmonics and their power spectra of the coefficients af m are composed of signal 
parts Cf p and noise parts Nf p : 



/ p * p' \ 

\ a lm a l'm'/ 



pp' 



(9) 



where PP' represents three pairs of maps, TT, EE and TE. The signal parts are computed 
from a fiducial cosmology. We assume the cosmological parameters of the WMAP5 mean 
values for a power-law ACDM model as a fiducial model. As for the noise power spectra, 
we assume a Gaussian beam and a spatially uniform Gaussian white noise. Nf p are given 
as the combined effects from these and can be approximated as 



pp' 



5 PP ,{d 



FWHMO- P ) 2 exp 



1(1 



, FWHM 

81n2 



(10) 



where #fwhm is the full width at half maximum of the Gaussian beam and a p is the 
root mean square of the instrumental noise. For the expected data from the Planck 
experiment, we use three frequency channels at 100, 143 and 217 GHz. We adopt the 
following values for the instrumental parameters [29]: (9 FW n W [aicm.m], cr T [/iK], <7p[/iK]) = 
(9.5, 6.8, 10.9), (7.1, 6.0, 11.4) and (5.0,13.1,26.7) for v = 100, 143 and 217GHz, respec- 
tively. We assume other frequency channels are used to remove foregrounds and they are 
ideally removed. 
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Since the anisotropies from both signal and noise are Gaussian distributed, the likeli- 
hood function of the data a = {af^} for a theoretical model with parameters = {^} 
is given by 

£(a|0) oc ^7= ex P (-"^[(7(0)]-^ , (11) 

where C(0) is a covariance matrix of the theoretical data. Denoting a covariance matrix 
of mock data as C, the effective \ 2 is given as 



Xc 2 ff = ]T(2/+l)/ sky 



Inigl + C-^-i 



(12) 



We take f s ^ y = 0.65 as the expected sky-coverage for the Planck experiment. The factor 
(21 + l)/sk y represents the effective number of independent moments obtained from the 
observation. For the MCMC analysis, we include the data up to / = 2500. 

Before presenting our results, here we comment on possible contributions from the 
thermal and kinetic SZ effect. We assume that the thermal SZ effect can be precisely 
estimated from the other lower frequency channels of the Planck survey than those used 
in our analysis, and can be removed ideally. The contribution from the kinetic SZ effect 
on CMB anisotropy depends on the details of the reionization process. For somewhat 
conventional scenario, as argued in [30] , it would be only about a few percent in the range 
of multipoles we make use of, I < 2500, and also sufficiently smaller than the expected 
instrumental noise for the Planck survey. Ref. [31] argued that "patchy" reionization would 
make it significantly larger but they found that the shape of the power spectrum due to the 
kinetic SZ effect does not depend much on the reionization model. Then, they concluded 
that its effect on the determination of the cosmological parameters can be neglected by 
marginalizing over the amplitude of the kinetic SZ power spectrum. We thus neglect the 
kinetic SZ effect here. 

Now we discuss a future constraint on Y p from the Planck experiment. In Fig. [31 one- 
dimensional marginalized likelihood for Y p is shown. For comparison, we also plot the 
constraint from current observations. We can expect that the uncertainty for Y p at 68% 
C.L. becomes as AY p ~ 10 -2 , which is 10 times smaller than that from current data (see 
Table [3] below). Since Planck can measure the CMB power spectrum at higher multipoles 
very precisely, the effects of damping due to Y p can be well probed. It should also be noted 
that since likelihood functions for Y p and other cosmological parameters have almost the 
Gaussian form, our results here using MCMC approach are almost the same as those 
obtained by Fisher matrix analysis which has already been done in Refs. [10,12]. Thus 
we found that the Fisher matrix analysis can give a good estimate for Planck data for 
the parameter set we assumed here. The results here are consistent with those given in 
Ref. [29,32] in which a forecast on Y p is investigated using MCMC analysis too. 

Next we discuss the effects of prior on Y p on the constraints on other cosmological 
parameters in the Planck experiment. As mentioned in the Introduction, when one tries 
to constrain some cosmological parameters from CMB, the value of Y p is fixed to be 
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Figure 3: One- dimensional marginalized distributions of Y p . Shown are the distributions 
from the Planck experiment for the cases with no priors on Y p (red solid line) and assuming 
the BBN relation (green dashed line). For reference, the distribution for the case with no 
priors on Y p using current CMB data is also shown (dash-dotted magenta line). 



0.24 in most of analysis. In the previous section, we showed that, when we use current 
cosmological data, the fixing of Y p = 0.24 does not affect much the constraints on other 
cosmological parameters since the value of Y p itself is not constrained well. However, as 
just shown above, Planck can measure the value of Y p much more precisely, thus we should 
study the effects of the assumption of Y p when we constrain other parameters. For this 
purpose, we made MCMC analyses for three cases: (i) Y p is not fixed but varied freely, (ii) 
Y p is fixed as Y p = 0.24 and (iii) Y p is regarded as a function of Ub via the standard BBN 
calculation. For the case (iii), we relate the value of Y p to cjj by the fitting formula given 
in [33], to which we refer as "BBN relation" in the following. 



parameters 


Y p free 


SBBN Y p {uj h ) 


Y = 0.24 


u b 
r 

Tig 

ln(10 10 A s ) 


U.U^Z (O_ 000 25 

1 nQ8 +a0015 
u.iuyo_ 0014 

1.040631SJ 

0.0879±8:gg__ 

OR97+0.0079 
U.JUZ 1 „ .0085 

3.064_ffl 


n n997q+ u - UUU17 

U.UZZ (O_ 00017 

inqq+ 0015 
u.iujy„ 0014 

1-04061+™ 

o.o88oi858g 

OfiQ 1 +0.0046 

u.yuoi_ 0042 
3.065_°;°1° 
0.248586_°;™ 


02261 +u - UUUiti 

o 1 1 no +0 - 0013 

u.iiuu„ 0016 

i - u ^ :UOi -0. 00038 

n n«7i +0.0049 

Q c i80 +0 - 0042 
u.yoou„ 00 4 4 

3 061 +a0l ° 

O.U01_o.oi2 


#0 


n 9^7+0.0080 
U.zou 1 _ .0086 

71.881°;*° 


u.zouj„ 0080 
71 Q9+ - 78 


n OKO7+0.UU74 
U.ZOCW _ .0083 

71 61 +0 ' 73 



Table 3: Mean values and 68% errors from Planck for some assumptions on Y p . 
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Figure 4: One-dimensional marginalized distributions of cj&, u c , n s , A s , r, H , using same 
data as in Fig. [3j Additionally, the distributions from Planck for the case with fixing 
Y p = 0.24 are also shown (dotted blue line). 



Now we show one-dimensional marginalized likelihood for Ub, u c , n s ,A s , r, H Q in Fig. HI 
In the figure, three cases (i), (ii) and (iii) are depicted. In table El the mean values and 
errors at 68% CL. are shown for representative parameters. By looking at Fig. HI some 
features can be noticed. For u c , A s and r, the effects of the prior on Y p are very small 
even with the precision of Planck. However, for u>b,n s and H , marginalized distributions 
are changed depending on the prior on Y p . This tendency can also be seen by reading the 
errors at 68% CL. from Table [31 For ujb,n s and H , when we assume the BBN relation 
or fix the value of the helium abundance as Y p = 0.24, the errors are reduced to some 
extent, which clearly indicates that the assumption of Y p can affect the determination 
of other cosmological parameters. Furthermore, for these parameters, when we fix Y p = 
0.24, the central values differ from the fiducial values by about the uncertainties at 68% 



C.Lcj Therefore, in the Planck era, we advocate varying the value of Y p freely in the 
cosmological parameter estimation for a conservative constraint, or, if we would like to do 
the cosmological parameter estimation in the framework of the standard cosmology, we 
should impose the BBN relation. 

To see how these parameters are correlated with Y p , 2D marginalized contours may 
be useful, which are shown in Fig. [51 From this figure, we can see that Y p and these 



^ 4 A similar analysis has been done in Ref. [32] recently and their results are consistent with ours. 
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parameters ujb,n s and Hq are positively correlated. Positive correlations of Y p with n s 
mainly come from the degeneracy at higher multipoles where the effect of the diffusion 
damping is significant. From Eqs. (13), (16]), (EJ) and (jHJ), the response of Hi for i > 2 to the 
change of Y p and n s are opposite sign, which indicates that the correlation between these 
are positive. Correspondingly, Ub and Y p becomes positively correlated because of the 
positive correlation between cj& and n s which can be read off, in particular from Eq. (j3J). 
For the correlation of Y p with H , it should be noticed that the position of the first peak 
can be significantly affected by changing n s and H . By increasing the value of Y p , the 
diffusion damping suppresses the power on small scales. To compensate this effect to fit 
the data well, increasing n s can enhance the power. Due to the change of n s , the first 
peak position is in turn also shifted. However, as is already well determined by WMAP, 
the position of the first peak should be tuned to the right position to fit the data well. 
This can be done by changing H as can be seen from Eq. (J3|). In fact, the change of Y p 
itself can also shift the peak position. However the direct effect of Y p on l\ is very small 
compared to other quantities. 




66 - 66 - 

0.022 0.023 0.024 0.022 0.023 0.024 0.94 0.96 0.98 1 

Figure 5: Two-dimensional marginalized constraints in several combinations of cosmo- 
logical parameters. Shown are a Planck forecast on the constraints for the cases with no 
prior on Y p (orange and yellow shaded region), assuming the BBN relation (solid red line) 
and fixing Y p = 0.24 (blue dotted line). For reference, the constraint for the case with no 
prior on Y p using current CMB data is also shown (dash-dotted magenta line). 

Finally, we would briefly discuss how our arguments above are affected by considering 
the uncertainties arising from the recombination process. Several authors have claimed 
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that the uncertainties in the theory of recombination make some effects on the CMB power 
spectra and the determination of the cosmological parameters from CMB [34,35]. Although 
some detailed studies of the recombination modeling have been done [36-41], more devel- 
opments of the recombination modeling may be needed [35] . Since the primordial helium 
abundance, which we are studying in this paper, directly affects the recombination history, 
it may be interesting to check how the uncertainties in the recombination modeling affect 
the determination of Y p and other cosmological parameters in the future Planck experi- 
ment. For this purpose, we treat two numerical parameters, the so-called fudge factors 
Fh and bn e used in recfast [35,42,43] as free parameters to represent the uncertainties 
of the recombination modeling. Fh is introduced to fit the hydrogen recombination rate 
in the three-level approximation to the result from multi-level calculations. bn e is a fitting 
parameter for the recombination rate of Hel. We repeated a MCMC analysis including Fh 
and bn e and marginalized over these parameters with the top hat priors given in Table [TJ 
which are very conservative ones. In Table H] we summarize the resultant constraints on 
the cosmological parameters from expected Planck data. We also show the probability 
distributions for several cosmological parameters in Fig. |6j Table H] shows that, even if we 
adopt a very conservative prior on the fudge factors, the uncertainties of the recombination 
modeling which is represented by Fh and bn e do not significantly affect the determination 
of cosmological parameters in the Planck era. (Errors for some parameters are changed 
at most by 10%.) Therefore, we can say at least that the theoretical uncertainties of the 
recombination process which are discussed recently do not affect our previous discussions. 
However, since our analysis is limited, we need more understanding of the recombination 
process and more detailed analysis in order to reduce systematic errors in the helium 
estimation from future CMB data. 





Y p free 


SBBN Y p (u b ) 


parameters 


F H , b He fixed 


Fh, b He free 


Fh, bne fixed 


F H , b He free 


o s 

T 

n s 

ln(10 10 yl s ) 


U.UZZ(O_ 000 25 

0.1098iS 

1.04063+0™ 
08 79+ - 0054 

U.UO ( y_ .0062 
Ofi97+0.0079 
U ' JUZI -0.0085 

3.0641S 
0.248lffi 


n 997Q+0. 00025 
U.ZZrfO_ .00024 

0.10981^ 
1.040661™ 

9628+ ' 0091 

U.JUZO_ 0081 
O.UU^_ 014 

n 940+0.013 

U.Z^J_ .012 


n 0997^+0-00017 
u ' uzz '°-0 .00017 

0.10991S 
l-04061+o™ 
0.08801S 
Qfi31+°. 0046 

U.yUOl_ .0042 

3.065+°;™ 
0.248586+8™ 


0.0022731--- 

0.109818:8^ 
i.04064i°:8S 
0.088018:88^ 

9626+ - 0046 

U.JOZO_o 0050 
flfi4+0.012 
O.UU^„o.012 
n O/|Cc;qq+0.000087 
U.248£)8<j„ .000058 


Ho 


n 9^7+0-0080 
U.zoui „ .0086 

71.881™ 


2964+ u ' UU8i 

U.ZOU^-o.0088 
71 OI +0-86 
' i - ai -0.82 


U.ZdUO„ .0080 

71 Q9+ - 78 
' L - JZ -0.66 


U.ZOUO„ .0081 

71.891°:^ 



Table 4: Comparison for the cases with and without uncertainties from the recombination 
process being considered. In the first two columns Y p is treated as a free parameter and 
in the latter two columns the SBBN relation is assumed. 
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Figure 6: Comparison for the cases with and without uncertainties from the recombination 
process being considered. Shown are the cases with Y p being treated as a free parameters 
and the fudge factors being fixed to the standard values (red solid line), Y p being treated 
as a free parameters and the fudge factors being treated as free parameters (blue dotted 
line), Y p from the SBBN relation and the fudge factors being fixed to the standard values 
(green dashed line) and Y p from the SBBN relation and the fudge factors being treated as 
free parameters (magenta dot-dashed line). 



5 Summary 

We studied the constraint on Y p and the effects of the priors for Y p on constraining other 
cosmological parameters using current CMB data from WMAP5, ACBAR, BOOMERANG 
and CBI, and also from the future Planck experiment. After briefly reviewing the effects 
of Y p on CMB, we studied current constraints on the primordial helium abundance. We 
obtained the current limit on Y p from WMAP5 alone as Y p < 0.44 at 95%C.L., which is 
improved to be Y p = 0.25^°^;^ at 68% (95%) C.L. by adding the data of ACBAR, 
BOOMERANG and CBI around the damping tail. We have also considered how the prior 
of Y p can affects the constraints on other cosmological parameters using currently available 
data. We found that, at the present precision level of CMB measurements, the prior on 
Y p has little effect for determinations of other cosmological parameters. 

We have also investigated the future constraint from the Planck experiment. By per- 
forming a MCMC analysis, we derived an expected error for the helium abundance from 
the future Planck experiment and found that it will be well measured with the accuracy 
of AY p ~ (68% C.L.) in the Planck experiment, which is 10 times smaller values 
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compared with current data. Furthermore, it may be interesting to notice that this preci- 
sion is comparable to that obtained by HII region observations. As for the effects of the 
prior on Y p on the determination of other cosmological parameters, we found that, with 
the precision of Planck, the assumption on Y p can affect the constraints on other cosmo- 
logical parameters such as Ub and n s . In this respect, the prior on Y p can be important 
for determining the other parameters. In addition, the constraint on Y p from CMB itself 
can be an independent test from other methods such as using HII region observations. 

In the near future, we can have more precise measurements of CMB. Such upcoming 
data would give us more precise information of the primordial helium abundance. At the 
same time, it is necessary to study the effects of the helium abundance more rigorously in 
order to extract information of other cosmological parameters. 



Note added: While we were finishing the present work, Ref. [32] appeared on the 
arXiv, which has some overlap with our analysis on the constraints from the future Planck 
experiment. 
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